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^ In a recent approach in modelling a radiating relativistic star undergoing gravitational 

, collapse the role of the Weyl stresses was emphasised. It is possible to generate a model 

which is physically reasonable by approximately solving the junction conditions at the 
boundary of the star. In this paper we demonstrate that it is possible to solve the 
Einstein field equations and the junction conditions exactly. This exact solution con- 
tains the Priedmann dust solution as a limiting case. We briefly consider the radiative 
; I ■ transfer within the framework of extended irreversible thermodynamics and show that 

' relaxational effects significantly alter the temperature profiles. 
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1. Introduction 

H ; 

I Since the pioneering work of Oppenheimer and Snyder^ the final state of a star 

undergoing gravitational collapse has occupied a significant standing in relativistic 
astrophysics. An interesting extension of this problem is to include heat flow to a 
spherically symmetric model as a realistic model of a radiating star in gravitational 
collapse. Although this model is simplified, it does provide a good insight into the 
nature of collapse and serves as an excellent yardstick for more complicated radi- 
ating models. The available literature^^-'^^ on radiative gravitational collapse, and 
relativistic astrophysical applications of spherically symmetric gravitational fields 
with an imperfect matter distribution, is extensive. Models of relativistic radiating 
stars are useful in the investigation of the Cosmic Censorship hypothesis^^"^^ which 
has attracted much attention in recent times amongst researchers in the field. 

Since the derivation of the junction conditions'* for a radiating spherically sym- 
metric star undergoing gravitational collapse and dissipating energy in the form of a 
radial heat flux, there have been many novel approaches to finding physically viable 
models. The model proposed by de Oliveira et al. ^ in which the star collapses from 
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an initial static configuration has been studied extensively. Recently^ it has been 
shown that a similar approach can be adopted to look at the end state of collapse 
which results in a superdense star. Another model which has been useful in under- 
standing the effects of dissipation is due to Kolassis'' in which the fluid trajectories 
are assumed to be geodesies. Recently Herrera et al ^ have imposed the condition of 
conformal flatness for a radiating, shcar-frcc mass distribution undergoing gravita- 
tional collapse. This approach generalises earlier works by Patel and Tikekar ^ and 
Grammenos^^ in particular. In their paper Herrera et al.^ provided an approximate 
solution of the Einstein field equations. Their approximate solution highlighted the 
relaxational effects of the collapsing fluid on the temperature profile. These results 
were in agreement with earlier work carried out on the thermodynamics of radiating 
stars within the framework^ of extended irreversible thermodynamics. In this 
paper wc demonstrate that it is possible to solve the field equations exactly with 
the assumption of vanishing Weyl stresses. Our exact solution provides a mathe- 
matical basis for the approximate approach of Herrera et al.^ and indicates that 
their model is physically signiflcant. 

This paper is organised as follows. In section two, we present the field equa- 
tions governing the interior spacetime comprising a shear-free, spherically symmet- 
ric matter distribution with a radial heat flux. We also impose the condition for 
conformal flatness. In section three, the junction conditions required for the smooth 
matching of the interior spacetime to the Vaidya spacetime are outlined. In section 
four, an exact solution to the boundary condition is presented. We show that our 
exact solution has a Friedmann limit as in the case of the approximate analysis of 
Herrera et al.^ when the heat flux is absent. The stability of the model is also briefly 
discussed. Section five deals with the evolution of the temperature profiles within 
the framework of extended irreversible thermodynamics. We are in a position to 
calculate the temperature explicitly in both the causal and noncausal theories and 
it is further shown that relaxational effects lead to higher temperature within the 
stellar core. 

2. Shear-free spacetimes 

We investigate the gravitational collapse of a shear-free matter distribution with 
spherical symmetry. This is a reasonable assumption when modelling a rclativistic 
radiating star. In this case there exists coordinates for which the line element may 
be expressed in a form that is simultaneously isotropic and comoving. With the 
coordinates (x") = (t, r, 9, (p) the line element for the interior spacetime of the 
stellar model takes the form 

ds^ = -A^dt" + [dr^ + r\de^ + sin^ Odcj)'')] . (1) 

where A = A{t, r) and B = B{t, r) are metric functions. In this paper we consider a 
model which represents a spherically symmetric, shear-free fluid conflguration with 
heat conduction. For our model the energy-momentum tensor for the stellar fluid 
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becomes 

Tab = {P + P)UaUb + PQab + (la^b + QbUa ■ 

The fluid four-velocity u is comoving and is given by 

The heat flow vector takes the form 

9"* = (0,3, 0,0). 



(2) 
(3) 
(4) 



since q°'Ua = and the heat is assumed to flow in the radial direction on physical 
grounds because of spherical symmetry. The fluid collapse rate = u?„ of the stellar 
model is given by 



e 



(5) 



where dots represent differentiation with respect to t. The nonzero components of 
the Einstein's fleld equations for the line element (1) and the energy-momentum 
(2) are: 

/2 
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9 = - 



AB"^ 



B' B'B A' B' 



(8) 



This is a system of coupled partial differential equations in the variables A,B,p,p 
and q. The nonvanishing components of the Weyl tensor Cabcd for the line element 
(1) are 



C2323 = jB^sin'^e 



A^ 
~A 



B^ 
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A^ 
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= 2r2(^^^ (7o303 = -2r2sin2^Ci2i2 = -2r2Ci3i3. (9) 

If we now impose the condition of conformal flatness, then this would imply van- 
ishing of all the Weyl tensor components. From the above relations we note that 
this condition is fulfilled if we demand 

A" B"\ [A' B'\ f B'\ 

Equation (10) can be easily integrated to give 

A{t,r) = [Ci{t)r'' + l]B (11) 

where Ci {t) is an arbitrary fiuiction of t yet to be determined. 
The condition of pressure isotropy reduces to 



which has the general solution 



^-C2(t)r2 + C3(i) ^^^^ 



where C2{t) and Cz{t) are arbitrary functions of integration. 



3. Junction Conditions 

Since the interior is radiating energy, the exterior spacetime is described by Vaidya's 
outgoing solution given by^^ 

= - (l - dv^ - 2dvdR + {de^ + sin2 Od^^) . (14) 

The quantity m{v) represents the Newtonian mass of the gravitating body as mea- 
sured by an observer at infinity. The metric (14) is the unique spherically symmetric 
solution of the Einstein field equations for radiation in the form of a null fluid. The 
Einstein tensor for the line element (14) is given by 

The energy-momentum tensor for null radiation assumes the form 

Tab = ^WaWb (16) 

where the null four-vector is given by Wa = (1, 0, 0, 0). Thus from (15) and (16) we 
have that 
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for the energy density of the null radiation. Since the star is radiating energy to the 

dm 

exterior spacetime we must have — — < 0. 

dv 

The necessary conditions for the smooth matching of the interior spacetime to 
the exterior spacetime was first presented by Santos^ in his seminal paper. The 
junction conditions for the line elements (1) and (14) are given by 

(rB)s = i?s (18) 
PS = (qBh (19) 

Matching the interior spacetime (1) to the outgoing Vaidya solution (14) leads to 
the following ordinary differential equation 

+2 (^^^) [C2(C2 - CiC3)62 + C3(CiC3 - 2C2)] = . (21) 

In the above Ci = Ci{t), C2 = C2{t) and C3 = C3{t). Dots denote differentiation 
with respect to t. Also rs = & is a constant which defines the boimdary of the star. 
To complete the gravitational description of this radiating star we must solve the 
junction condition (21). 



4. An exact model 

Herrera et al.^ have presented a simple approximate solution of (21) which reduces 
to the Friedmann dust sphere in the appropriate limit. They assumed the following 
forms for the temporal functions 

Ci=eCi{t), C2 = 0, ^3 = ^. (22) 

where < e << 1 and a is a positive constant. With these assumptions, the junction 
condition (21) reduces to 

C^i+(^ + ^)^?i-0 (23) 
which has the approximate solution 

Ci«Ci(0)exp(^-|). (24) 

The approximate solution (24) makes it possible to perform a qualitative analysis 
of the physical features of the model. 
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However it is possible to solve (21) exactly for particular choices of the arbi- 
trary temporal functions. Wc now present a simple exact solution of the junction 
condition (21). Guided by the approximate solution of Herrera et al.^ we take 



C2(i) =0. 

Substituting (25) into the boundary condition (21) we obtain 



Ci62 + 1 



2C\b 



C3C3 + 2{CW + l)CiCi = . 



The transformation 

enables us to write (26) in the equivalent form 



u — u 



Clip + 1 



+ 2C16 



-(Ci6^ + l)Ciu = 0. 



(25) 



(26) 



(27) 



(28) 



Equation (28) is a second order differential equation which is linear in u{t) if the 
function Ci is specified. A simple choice for Ci is 

Ci=C (29) 

where C is a constant. It follows that (28) becomes 

u - 2Cbu - {Cb^ + l)Cu = . (30) 

This differential equation has three categories of solutions. 

In terms of the original function C3 we can present the solutions as 
Case I: 20^"^ + C > 



Case II: 2C%'^ +C<0 



(31) 



Csit) = /Jie'^" cos {tV2C%'^+C) + p2e^''* sin {tV2C%^ + C) 
Case III: 2C%^ +C = 



(32) 



C3{t) = {Pi+l32t)-'e 



2„-2Cbt 



(33) 



where (3i and /?2 are constants of integration. We have therefore exhibited an exact 
solution to the boundary condition (21) where 



Ci{t)=C, C2{t) = 0, 



(34) 
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and the three analytic forms for C3 are given by (31)-(33). This exact solution 
provides a formal mathematical basis for the approximate model of Herrera et al.^ 
for a radiating star. 

We now consider Case III more closely. The line element in this case becomes 



ds^ = (/3i + P2t)*e*'''" I - (Cr^ + ifdf + dr^ + r^d0^ + sin^ 0d<P^)\ . (35) 

If we set /3i = C = 0, /J2 = 1 in (35) we obtain 

ds^ = [-dt^ + dr^ + (de^ + sin^ ed<l)^)] (36) 

which is the Friedmann dust solution. Herrera et al.^ also regain (36) for their 
approximate solution. However our limit arises from an exact model. 
From (6)- (8) we obtain the expressions 



12 [/32 + bCzf e-4CW 



(Cr2 + 1)2 



(37) 



4C [(1 - h'^C + Cr^)pi + {-2b + t+{r^- b^)Ct)p2] 

P = ^ WTW ^''^ 

for the line element (35) where z{t) = /3i + (52t- The luminosity radius is given by 

(rS)i: = (/?i+/?2i)'e2^''* (40) 

which starts off at arbitrary large values at t = —00 and evolves towards t = Q. The 
four-acceleration of this model is nonvanishing as 

A' ICr 

a. = .''V,.. = -61 = (41) 

which is also independent of time. The measure of the dynamical instability of the 
stellar configuration at any given instant in time is given by 

r = ^. (42) 

alnp 

For our model we obtain 

r -r _ h^c\i3, + [i2tf 

^ centre ^ surface - 6(/32 + 6C(/3i + /32i)' ^ ^ ' 

This implies that Tcenter < ^surface for all times. This physical result indicates that 
the centre of the collapsing star is more unstable than the outer regions, reinforcing 
the approximate results of Herrera et al.^ as well as the earlier work by Denmat et 
alP^ which is also related. 
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5. Thermodynamics 

In this section we investigate the evolution of the temperature profile of our model 

within the context of extended irreversible thermodynamics. The causal transport 
equation in the absence of rotation and viscous stress is 

Tha''qb + qa = -K {ha^VbT + Tiia) (44) 

where hab = gab + UaUb projects into the comoving rest space, T is the local equi- 
librium temperature, k (> 0) is the thermal conductivity, and r (> 0) is the rclax- 
ational time-scale which gives rise to the causal and stable behaviour of the theory. 
To obtain the noncausal Fourier heat transport equation we set r = in (44). For 
the metric (1), equation (44) becomes 

r{qB) + AqB = -'^^^ . (45) 

In order to obtain a physically reasonable stellar model we will adopt the thermo- 
dynamic coefficients for radiative transfer. Hence we are considering the situation 
where energy is transported away from the stellar interior by massless particles, 

moving with long mean free path through matter that is effectively in hydrody- 
namic equilibrium, and that is dynamically dominant. Govender et af'^ have shown 
that the choice 

is a physically reasonable choice for the thermal conductivity k, the mean collision 
time between massive and massless particles Tc and the relaxation time r. The 
quantities a > 0, P > and a > are constants. Note that the mean collision time 
decreases with growing temperature as expected except for the special case u = 0, 
when it is constant. With these assumptions the causal heat transport equation 
(45) becomes 

rp3-<7( Aj^y 

(i{qB)T-^ + A{qB) = -a . (47) 

Solutions are easily obtainable when in the noncausal case (/3 = 0) as was shown 
by Govinder and Govender^^ in their general treatment. All noncausal solutions of 
(47) are 

{Aff-'' = / A^-'^qB^dr + Fit), a (48) 

a J 

In {AT) = -- I qB'^dr + F{t), a = 4 . (49) 

a .1 



where F{t) is a function of integration which is fixed by the surface temperature of 
the star. Note that T corresponds to the noncausal temperature when /3 = 0. For 
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a constant mean collision time (cr = 0), (47) can be integrated to give the causal 
temperature, ie., 

4 " 



(AT)^ = -- 



J A^B{qB)dr + j A^qB'^dr 



(50) 



In (46)we can think of /3 as the 'causality' index, measuring the strength of relax- 
ational effects, with (3 = giving the noncausal case. 
The effective surface temperature of a star is given by 

1 \f±_\ 



(51) 



^,2^2 J y^T^S J 

where L is the luminosity at infinity and 6{> 0) is a constant. The luminosity at 
infinity can be calculated from 



dm 

J-'oo = — r' 
dv 



where 



m{v) 



2^2 



3r/2 



r'B' 



r-^B 



2B 



(52) 
(53) 



For our model we can calculate the temperature in both the causal and noncausal 
theories explicitly. The noncausal temperature for constant collision time is given 

by 



1 + C62 



1 + Cr2 



(1+Cr2)3-(1+ C6^)3 



(54) 



where z{t) = /3i + /32i- The causal temperature for constant collision time is given 
by 



rp4 



'1 + Cb^' 






_1 + Cr2_ 


a 





+ 



16 {p2Z-'^ + Cb) 
3a 



z 



(1+Cr2)2- (1+^2)2 



(1 + Cr2)'^ 



(1 + Cr2)^-(1 + C62) 



,2\3 



-2C6t 



(55) 



(1 + Cr2)4 

We note that at the boundary of the star, the causal and noncausal temperatures, T 
and T respectively, are equal. In fact the causal temperature is everywhere greater 
than its noncausal counterpart within the interior of the star. From (45) we can 
write 



k{T){AT)' - K{f){Afy = -{BT){qB) 
which takes the form 



(56) 
(57) 



for the assumptions (46) where we have defined T = AT and T = AT. We can 
immediately see that the relative spatial gradient of T is greater than that of T, 
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since {qB) > and 4 — a > 0. This result confirms earlier findings^^ in which the 
acceleration was vanishing and the model had a Friedmann limit. 

The evolution of the temperature has been studied within the framework of 
extended irreversible thermodynamics and are in keeping with earlier works on 
radiating collapse. This model has shown that relaxational effects can significantly 
alter the temperature of the star, especially during the latter stages of collapse. We 
have further shown that the model is more unstable closer to the centre as opposed 
to the outer regions of the star. Such instablities have been shown to result in 
peeling or cracking up^^ of the stellar configuration. 
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